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ABSTRACT 

Scour  About  A  Single,  Cylindrical 
Pile  Due  to  Random  Waves  and  Current.  (December  1984) 
Robert  W.  Eadie  IV,  B.S.,  United  States  Naval  Academy 
Chairman  of  Advisory  Committee:   Dr.  J.  B.  Herbich 

An  empirical  approach  was  utilized  in  this  investiga- 
tion to  study  the  scouring  effects  that  various  currents 
and  random  waves  have  on  a  single  cylindrical  pile.   Re- 
sults from  laboratory  experiments  are  presented  in  graphi- 
cal form  depicting  relationships  between  pertinent  dimen- 
sionless  parameters  which  were  developed  through  a  dimen- 
sional analysis.   These  results  were  then  compared  to  re- 
sults compiled  in  previous  studies  that  investigated  scour 
due  to  currents  and  monochromatic  waves  in  order  to  deter- 
mine if  the  irregularity  of  the  wave  train  is  significant 
in  the  study  of  scour. 
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PART  I 
INTRODUCTION 

The  construction  of  offshore  structures  is  an  ex- 
tremely expensive  process.   A  portion  of  these  construc- 
tion costs  are  due  to  large  safety  factors  used  during 
the  design  process.   These  safety  factors  which  may  lead 
to  overdesign  and,  therefore,  unnecessary  expense  are  the 
result  of  the  engineer's  inability  to  accurately  predict 
some  of  the  forces  on  the  structure,  as  well  as  some  of 
the  natural  phenomena  that  will  occur  once  the  structure 
is  in  place.   One  of  these  natural  phenomena  which  is  be- 
coming an  ever  increasing  problem  is  scour  caused  by  waves 
and  currents  around  offshore  structures.   Design  factors 
for  offshore  structures  depend  on  the  location  of  the  ocean 
bottom.   The  scour  of  a  non-cohesive  ocean  floor  can  lower 
the  bottom  line  significantly,  and  by  doing  so  can  dras- 
tically change  these  design  factors. 

The  purpose  of  this  investigation  was  to  study  what 
effect  various  currents  and  random  waves  have  on  a  singu- 
lar, cylindrical  pile.   These  results  were  then  compared 
to  results  compiled  in  previous  studies  that  investigated 
scour  due  to  currents  and  monochromatic  waves  to  determine 
if  the  irregularity  of  the  waves  is  significant.   An  em- 
pirical approach  was  used  in  this  study. 


Experiments  were  conducted  in  a  two-dimensional  wave 
tank.   The  pile  used  in  this  study  had  a  diameter  of  1.5 
inches.   Two  water  depths,  two  currents,  one  sediment  size 
and  four  random  wave  spectra  were  utilized.   Using  data 
obtained  from  the  experiments,  an  attempt  was  made  to  de- 
scribe scour  in  terms  of  relevant  dimensionless  and  dimen- 
sional parameters.   By  doing  so,  it  is  intended  that  the 
reader  will  obtain  a  better  understanding  of  the  scouring 
process . 


PART  II 
REVIEW  OF  LITERATURE 
Incipient  Motion 

The  scour  phenomenon  begins  with  the  incipient  motion 
of  sediment  particles.   Eagleson  and  Dean  (7)  define  in- 
cipient sediment  motion  to  be  "an  instantaneous  condition 
reached  when  the  resultant  of  all  the  active  forces  on 
the  particle  intersects  the  line  connecting  the  bed  parti- 
cle contact  points."   There  are  three  types  of  active 
forces  acting  on  a  sediment  particle.   Two  of  these  forces, 
the  drag  and  lift  forces,  are  due  to  the  fluid  motion. 
These  forces  are  in  turn  balanced  by  the  sediment  parti- 
cle mass.   These  forces  are  depicted  in  Figure  1  (1). 

A  thorough  understanding  of  these  forces  is  required 
in  order  to  comprehend  fully  the  processes  of  incipient 
motion,  bed-load  transport  and  bed  form  development;  all 
which  depend  on  the  collective  response  of  initially  sta- 
tionary bed  particles  to  fluid  forces.   The  following  dis- 
cussion of  incipient  motion  is  for  a  unidirectional  flow. 

The  variation  of  drag  forces  with  flow  condition  has 
been  well  documented,  as  has  the  lift  force  in  hydrodyna- 
mically-rough  flow  conditions.   However,  investigations 
of  lift  forces  in  hydrodynamically  smooth  and  transitional 
flow  have  yielded  conflicting  results. 
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Fig,  1.  -  Hydrodynamic  force  on  an  individual  sediment 
particle  (1  ) 


Prior  to  1967,  most  investigations  measuring  fluid 
lift  forces  were  conducted  in  the  hydrodynamic  rough  flow 
region  characterized  by: 

R*  =  V^D   >  70 


in  which  R*  represents  the  grain  Reynolds  number;  VVc  = 
shear  velocity;  D  =  grain  diameter;  and  v  =  kinematic  vis- 
cosity of  the  fluid.   All  measurements  under  these  hydro- 
dynamic  conditions  indicated  positive  lift. 

Coleman  (4)  was  the  first  to  study  lift  forces  in 
the  hydrodynamically  smooth  region  (R^  <  5)  and  transition- 
al flows  (5  <  R^.  <  70)  .   His  work  reported  measured  drag 
forces  and  inferred  lift  forces  from  observing  the  inci- 
pient motion  of  spheres  having  various  densities.   From 
the  results,  it  was  suggested  that  for  a  sphere  resting  on 
top  of  a  closely  packed  bed  of  spheres  of  equal  diameter, 
the  lift  force  decreased  through  zero  and  became  negative 
(bedward-directed)  as  Rv.  decreased  through  a  value  of  ap- 
proximately 15. 

Davies  (6)  prepared  an  experiment  whose  sole  purpose 
was  to  define  the  flow  conditions  under  which  negative 
lift  forces  occur.   He  found  that  significant  negative 
lift  can  be  expected  to  occur  only  if  the  following  con- 
ditions are  met:   1)  The  value  of  R,v  <  5  and  2)  The  test 
sphere  is  situated  such  that  a  significant  volume  rate  of 
flow  can  occur  beneath  the  sphere.   In  a  situation  in 


which  considerable  fluctuation  of  bed  shear  stress  occurs, 
the  onset  of  significant  negative  lift  with  a  decrease  in 
R*,  will  occur  gradually  and  unsteadily  as  the  effect  of 
turbulence  in  the  area  of  the  bed  particle  becomes  less 
significant,  until  hydrodynamically  smooth  flow  is  es- 
tablished at  Rj.  <  5  and  the  lift  force  becomes  predominant- 
ly negative. 

The  existence  of  varying  types  of  lift  forces  is  im- 
portant with  regard  to  bedform  development.   If  the  par- 
ticle is  removed  from  its  position  of  rest  at  R.v  <  5,  the 
negative  lift  forces  acting  on  the  particle  will  tend  to 
prevent  the  particle  from  leaving  the  immediate  vicinity 
of  the  bed,  and,  therefore,  it  will  move  along  the  bed  as 
a  bed  load.   On  the  other  hand,  if  the  particle  is  removed 
at  L  »  5,  it  will  experience  a  positive  lift  force  that 
will  tend  to  impart  an  upward  acceleration  and  thereby 
lift  the  particle  into  the  main  turbulent  flow,  where  it 
may  or  may  not  become  suspended. 

There  is  considerable  evidence  to  suggest  that  a  rip- 
pled bedform  exists  only  at  L  <  5,  and  that  a  dune  bedform 
exists  on  its  own  only  at  R*  >  24.   It  is  therefore  sug- 
gested that:   1)  The  flow  conditions  that  are  necessary 
for  existence  of  the  negative  lift  forces  on  a  sediment 
particle  correspond  closely  to  those  required  for  a  plane 
or  rippled  bedform.   2)  Positive  lift  force  conditions 
correspond  to  a  dune  bedform. 


The  condition  of  incipient  motion  for  a  cohesionless 
soil  particle  is  described  in  terms  of  the  forces  acting 
on  the  particle  by: 

F1 

F 


«  TT- (1) 


n 

where  F   and  F  are  the  forces  parallel  and  normal  to  the 

angle  of  repose  0.   In  our  case,  F.  and  F   are  the  resultant 
°  r  .        '   t      n 

forces  of  the  three  active  forces  mentioned  earlier,  the 
hydrodynamic  drag  force  F,,  the  lift  force  F,  ,  and  the 
submerged  weight  W  (gravity  force) .   In  accordance  with 
equation  (1) ,  the  condition  of  incipient  motion  under  the 
influence  of  these  three  forces  becomes: 

Wsina  +  F, 

tan  0  =  n ir=- (2) 

Wcosa  -  F, 

where  the  angle  a  is  the  inclination  of  the  bed  from  the 
horizontal  at  which  the  incipient  motion  occurs.   The  above 
situation  is  illustrated  in  Figure  2. 
tt  2 

fd  -  SKid     -r- (3) 

and 

tt   2 

fl  ■  cLK2d     -r- (4) 


where 


U,  =  bottom  fluid  velocity 

D 

Cn  =  drag  coefficient 
C,  =  lift  coefficient 


g> 


Fig.  2,  -  Force  diagram  on  particles  in  a  cohesionless 
loose  bed  (9) 


d  =  particle  diameter 

K1,K'2  =  s^aPe  factor 
p  =  fluid  density 
The  submerged  unit  weight  of  the  particle  is  expressed 
as  : 

W  =  K3(ps  -  P)gd3 (5) 

where 

K~  =  shape  factor 

p   =  particle  density 
By  introducing  equations  (3) ,  (4)  and  (5)  into  equation 
(2),  the  result  becomes, 

2 
(U,  )  cr      2K-,  (tan  0  cos  a  -  sin  a) 

(ps/p  -l)gd  =   CDK1  +  CLK2  tan  0    (6) 

where  (U,  )cr  is  the  critical  bottom  velocity  at  which  ac- 
cording to  equation  (1)  incipient  motion  occurs. 

(Ub)cr  =VV^ (7) 

where 

t   =  critical  shear  stress 
o 

p  =  fluid  density 
The  value  of  the  right  hand  side  of  equation  (6)  de- 
pends on  A)  The  particles:   their  size,  shape,  size  dis- 
tribution, uniformity,  .  .  . ;  B)  The  flow  dynamics,  which 
determine  the  values  of  CD  and  CL;  C)  the  channel  slope 
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which  in  many  cases  amounts  to  cos  a  ~  1;  and  D)  The  angle 
of  repose  which  depends  on  the  particle  properties  (9) . 

Once  the  critical  flow  velocity  is  reached,  the  cri- 
tical shear  stress  will  also  be  obtained  as  per  equation 
(7) .   Most  data  on  critical  shear  stress  for  non-cohesive 
sediments  have  been  observed  in  flume  experiments.   These 
experiments  have  shown  that  the  motion  of  the  sediment 
particles  at  the  bed  for  uni-directional  flow  is  highly 
unsteady  and  nonuniformly  distributed  over  the  entire  bed 
area.   Near  critical  conditions,  the  particle  motion  in 
any  small  area  of  the  bed  occurs  in  gusts  whose  incidence 
increases  as  the  mean  shear  stress  increases.   When  ob- 
serving a  large  section  of  a  sediment  bed  as  the  shear 
stress  approaches  the  critical  value,  the  occurrance  of 
these  gusts  of  particle  motion  seem  to  be  random  in  both 
time  and  space. 

Two  of  the  best  known  and  widely  accepted  investiga- 
tions into  the  concept  of  critical  shear  stress  are  those 
of  Shields  and  White  (17).   Shields  developed  the  relation: 

t          U^d 
7 2—^  =  f  (— ) (8) 

where 

t   =  critical  shear  stress 
o 

Y  =  specific  weight  of  sediment 


y  =  specific  weight  of  fluid 
d  =  particle  diameter 
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U.v  =  shear  velocity 

v  =  kinematic  viscosity 

Shields'  results,  in  which  he  determined  the  function 

U*d 
f  ( )  is  shown  in  Figure  3. 

As  discussed  earlier  under  Davies  experiment,  under 
flow  conditions  in  which  R^  <  5,  the  particle  will  have 
negative  lift  and  will  therefore  be  less  likely  to  move. 
This  corresponds  to  the  Shields  curve  that  indicates  that 
when  Rj.  decreases  below  about  10,  particles  increasingly 
resist  the  dimensionless  shearing  stress  of  a  flow. 

White  (17)  determined  the  critical  shear  stress  for 
several  sediments  under  conditions  of  developed  laminar 
flow,  laminar  boundary  layer  flow  and  turbulent  boundary 
layer  flow.   White  concluded  that  the  true  critical  shear 
stress  necessary  to  move  a  sediment  particle  has  a  fixed 
value  and  that  this  value  is  given  by  experiments  dealing 
with  laminar  flow.   White's  (24)  equation  for  x   for  sedi- 
ments in  a  horizontal  bed  is 

t   =  0.18  (v   -  y)d     tan  0 (9) 

c  s    '   s 

where 

2 
x   =  critical  shear  stress  (dynes/cm  ) 
c 

0  =  the  angle  of  repose  of  the  sediment  submerged 
in  the  fluid 
and  in  which  the  constant  (0.18)  is  obtained  from  experi- 
ments with  laminar  flow. 
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Fig,   3»   -   Shields   diagram  with  White   data  added    (2k) 
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Yang  (28)  expressed  some  reservations  on  using  Shields' 
diagram  as  the  criterion  for  incipient  motion.   Using  basic 
concepts  in  fluid  mechanics,  the  boundary  layer  theory 
and  laboratory  and  field  data  from  various  investigators, 
Yang  plotted  the  relationship  between  dimensionless  criti- 
cal velocity  and  shear  velocity  Reynolds  number  at  inci- 
pient motion  (Figure  4).   From  this  curve,  the  following 
equations  of  incipient  motion  were  obtained: 

vrT?        ?  5  U*d 

CR  =  £J> +  o.66;  for  0  <  — -  <  70  .  .  .   (10) 

log  (-2L-)  -0.06 


2.05;  for  70  <  (11) 


_CR 

W 

where 

V   =  average  critical  velocity  at  incipient  motion 
CR 

W  =  terminal  fall  velocity,  which  is  reached  when 
the  particle  drag  force  and  the  submerged  unit 
weight  of  the  particle  are  balanced. 

d  =  particle  diameter 

v  =  kinematic  viscosity 
U  =  shear  velocity  (IL.  =  ^ I    ,   P  =  fluid  density, 
t  =  average  shear  stress) 

Uu.d 
~   =  shear  velocity  Reynolds  number 

v 

Many  experiments  have  shown  that  for  any  given  bed 
material,  there  is  no  single  value  of  bed  shear  stress  or 
velocity,  below  which  not  a  single  particle  will  move  and 
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pig.  /+.  -  Relationship  between  dimensionless  critical 
velocity  and  shear  velocity  Reynolds  number 
at  incipient  motion  (23) 
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above  which  all  the  particles  having  the  same  size  will 
move.   The  forces  acting  on  the  particle  are  not  constant, 
but  instead  fluctuate  about  some  average  value,  a  phenomenon 
closely  associated  with  turbulence.   Also,  the  forces  which 
move  the  particle  are  dependent  on  the  particle  position. 
Therefore,  as  mentioned  earlier,  the  movement  of  bed  ma- 
terial particles  is  very  unsteady,  and  non-uniformly  dis- 
tributed over  the  entire  bed  (17)  .   Einstein  (8)  has  sug- 
gested that  the  process  of  bed  material  movement  is  sta- 
tistical in  nature.   Because  of  this,  it  appears  that  no 
truly  critical  condition  for  particle  movement  exists. 
This  is  further  exemplified  in  Figures  5  and  6  which  show 
the  variance  in  values  for  the  critical  shear  stress  and 
critical  water  velocity  for  quartz  sediment. 

Scour 

Fluid  Mechanics  of  Scour 

Local  scour  is  defined  as  "the  abrupt  decrease  in 
bed  elevation  near  a  pier  due  to  erosion  of  bed  material 
by  the  local  flow  structure  induced  by  the  pier."  (22) 
One  type  of  classification  of  the  type  of  scour  is  the 
amount  of  sediment  transport  into  and  out  of  the  scour  hole. 

Laursen  (12)  suggested  that  the  scour  process  could 
be  represented  by  the  relationship: 

^=q  ,  -  ,  , (33) 

3F^   4sl   Hs2 
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Fig.  5.  -  Critical  shear  stress  for  quartz  sediment  in 
water  as  a  function  of  grain  size  (2.1+) 
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Fig,  6.  -  Critical  water  velocities  for  quartz  sediment  as 
a  function  of  mean  grain  size  (2.1+) 
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in  which  dQg/dt  =  the  rate  of  local  scour  in  volume  per 
unit  time;  q  ^  =  the  capacity  of  the  flow  to  transport 
sediment  out  of  a  scour  hole  in  volume  per  unit  time;  and 
qs2  =  the  rate  at  which  the  undisturbed  flow  supplies  sedi- 
ment to  the  scour  hole. 

From  this  expression,  three  cases  of  scour  may  occur. 
They  include:   1)  no  scour,  in  which  case  q  ,  =  q  ?  =  0; 

2)  clear  water  scour,  where  q  ,  >  0  and  0  ~  q  9  «  q  , ;  and 

3)  scour  with  continuous  sediment  motion,  where  q  ,  >  q  0  >  0 

nsl    ns2 

Abrupt  changes  in  bed  elevation  near  a  pile  occur  only  in 
cases  (2)  and  (3) .  These  two  cases  of  scour  are  depicted 
in  Figure  7. 

The  main  flow  feature  responsible  for  scour  whenever  a 
slowly-varying  boundary  layer  current  is  present,  results 
from  secondary  flow.   The  passage  of  a  steady  flow  boundary 
layer  along  the  surface  past  the  cylinder  causes  a  stagna- 
tion pressure  gradient.   Because  the  stagnation  gradient 
decreases  "toward  the  sediment  bed,  a  secondary  flow  is 
generated  in  the  direction  of  the  bed.   This  secondary  flow 
has  two  effects  on  the  shear  or  boundary  layer  flow.   The 
principle  effect  is  the  formation  of  a  horseshoe  vortex 
that  wraps  itself  about  the  cylinder  and  trails  off  down- 
stream.  The  horseshoe  vortex  shown  in  Figure  8  is  depicted 
as  a  spiral  curve  extending  around  the  cylinder  base  (15) . 

Considerable  laboratory  data  are  available  from  scale 
model  studies;  however,  most  of  these  studies  investigate 
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Fig.  7.  -  Ti;ne  variation  of  scour  depth  (2.2.) 
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Fig.  3.  -  Secondary  flow  about  a  cylinder  (15) 
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the  horseshoe  vortex  in  a  river  environment.   Even  though 
these  results  provide  some  insight  into  the  characteris- 
tics of  flow  and  the  basic  mechanism  of  generation  of  the 
horseshoe  vortex,  it  is  impossible  to  apply  these  river 
model  results  directly  to  an  ocean  environment.   A  major 
deficiency  being  the  lack  of  appropriate  dynamic  and  kine- 
matic scaling  laws.   Shen  et  al .  (21)  have  shown  that  the 
size  and  the  strength  of  the  horseshoe  vortex  can  be  re- 
lated to  a  Reynolds  number  using  the  characteristic  length 
of  the  cylinder  in  the  flow  field.   This  relationship  is 
as  follows: 

a     2aU 
Vortex  strength  -  ££  -  f( -) (34) 

where 

co  =  rotational  velocity  of  vortex  core 
A  =  area  of  vortex  core 
a  =  half  width  of  pile 
U  =  undisturbed  flow  used  in  the  potential  flow  sense 

CO 

v  =  kinematic  viscosity  of  fluid 

Also,  Shen  et  al .  (21)  suggested  that  for  cylinders 
with  a  height  greater  than  the  thickness  of  the  bottom 
boundary  layer,  the  area  affected  by  the  horseshoe  vortex 
extends  approximately  one  cylinder  diameter  upstream  of 
the  obstruction. 

Hjorth  (10)  indicated  that  the  fluid  shear  stress 
striking  the  bottom  below  the  horseshoe  vortex  system  are 


increased  by  a  factor  of  eleven  when  compared  with  the 
shear  stress  values  associated  with  the  undisturbed  free- 
stream.   This  increase  in  bottom  shear  stress  causes  the 
bottom  sediments  to  scour.   With  scour  occurring,  the  bot- 
tom elevation  is  altered,  which  in  turn  brings  about  a 
change  in  the  flow  in  the  horseshoe  vortex. 

Another  effect  of  the  secondary  flow  about  a  cylindri- 
cal obstruction  is  the  downward  deflection  of  the  external 
flow  streamlines  about  the  cylinder.   This  downward  deflec- 
tion of  streamlines  is  shown  on  Figure  8,  the  secondary 
flow  is  depicted  at  the  stagnation  line  on  this  same  fig- 
ure.  The  convective  effect  of  the  streamlines  has  little 
effect  on  the  horseshoe  vortex,  except  for  a  slight  strength- 
ening of  the  rotational  vortex. 

A  circular  cylinder  blocks  the  approach  flow,  result- 
ing in  convective  acceleration.   The  acceleration  causes 
an  increase  in  the  bottom  shear  stress  by  a  factor  of  up 
to  four  over  the  bottom  shear  stress  of  an  unobstructed 
flow. 

Flow  patterns  in  the  cylinder  wake  also  influence 
scour.   There  is  a  net  downward  flow  in  the  wake  of  a 
cylinder  resting  on  a  plane  surface.   This  is  due  to  the 
fact  that  the  wake  region  is  a  low  pressure  region.   The 
wake  in  the  high  velocity  region  is  at  a  higher  pressure 
than  the  wake  in  the  low  velocity  region  near  the  sedi- 
ment bed.   This  difference  in  wake  pressure  causes  a  flow 
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to  form  in  the  wake  acting  downwards  toward  the  bed.   If 
the  effects  of  the  horseshoe  vortex  are  included,  this 
downward  flow  is  enhanced  (15) . 

In  an  attempt  to  get  a  better  understanding  of  the 
mechanisms  of  scour,  the  subject  will  be  dealt  with  in 
three  parts.   First,  the  mechanics  of  scour  in  a  steady 
current  alone  will  be  discussed,  followed  by  the  mechanics 
of  scour  due  to  oscillatory  flow,  and,  lastly,  the  mechanics 
of  scour  caused  by  the  interaction  of  waves  and  currents. 
Mechanics  of  Scour  in  Steady  Currents 

Repeated  observations  indicate  that  the  scale  and 
shape  of  the  scour  hole  closely  resemble  the  scale  and 
shape  of  the  horseshoe  vortex.   The  scour  is  caused  by 
the  large  bottom  fluid  shear  stress  gradients  occurring 
beneath  the  horseshoe  vortex.   The  initial  vortex,  before 
the  onset  of  scour,  exerts  bottom  fluid  shear  stresses 
that  are  approximately  one  order  of  magnitude  greater  than 
the  shear  stresses  in  the  surrounding  undisturbed  flow. 
Therefore,  it  is  possible  for  scour  to  occur  even  though 
sediment  is  not  eroded  and  transported  into  surrounding 
regions.   If  general  bedload  transport  is  occurring,  then 
the  elevated  bottom  shear  stresses  beneath  the  horseshoe 
vortex  produce  more  rapid  sediment  transport  and  bed  ero- 
sion. 

As  the  erosion  of  sediment  continues,  the  intensity 
of  the  horseshoe  vortex  decreases  as  it  spreads  out  into 
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the  scour  cavity.   As  the  vortex  intensity  decreases,  the 
fluid  shear  stress  acting  on  the  sediments  in  the  bottom 
of  the  scour  cavity  are  also  decreased.   Once  equilibrium 
is  reached,  the  fluid  shear  stress  acting  upstream  at  the 
bottom  of  the  scour  cavity  is  equal  to  the  fluid  shear 
stress  acting  downstream  on  the  undisturbed  portion  of 
the  bed.   The  fluid  shear  stress  acting  on  the  upstream 
face  of  the  scour  hole  is  equal  to  the  fluid  shear  stress 
acting  on  the  undisturbed  area  of  the  sediment  bed,  plus 
an  increment  approximately  equal  to  the  fluid  shear  stress 
required  to  entrain  the  sediment.   Once  this  distribution 
of  shear  stresses  has  been  achieved,  the  scour  hole  has 
reached  its  equilibrium  size  and  shape  (Figure  9)  and  the 
rate  of  transport  of  sediment  from  the  scour  hole  equals 
the  rate  of  supply  (15) . 
Mechanics  of  Scour  for  Oscillating  Flow 

Observations  of  scour  development  from  oscillating 
flow  are  given  by  Ranee  (19)  and  Palmer  (18) .   The  exact 
sequence  of  scour  hole  development  depends  on  the  diameter 
of  the  cylinder,  the  intensity  of  the  oscillating  flow, 
and  the  sediment  transport  on  the  surrounding  bed. 

Ranee  states  that  for  a  slender  body  (i.e.,  the  dia- 
meter of  the  body  is  small  when  compared  with  the  wave- 
length) under  the  action  of  waves,  the  shape  of  the  scour 
hole  is  dependent  on  the  amplitude  of  the  wave  motion  at 
the  bed.   When  the  amplitude  is  very  small  when  compared 
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Fig.    9.    -  Steady   current   ccour  cavity    (15) 
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to  pile  diameter,  two  discrete  scour  holes  are  formed  on 
the  forward  part  of  the  cylinder  at  45°  to  each  side  of 
the  center line.   For  this  case,  eddy  shedding  has  no  ef- 
fect on  sediment  transport.   When  the  bed  semi-orbit  length 
is  larger  and  becomes  approximately  one  fifth  of  the  pile 
diameter,  there  is  a  tendency  for  flow  separation  to  occur. 
Accretion  and  erosion  are  localized  and  four  scour  holes 
versus  two  are  formed.   These  scour  holes  have  a  maximum 
depth  of  about  10%  of  the  pile  diameter.   At  larger  ratios 
of  semi-orbit  length  divided  by  pile  diameter,  of  about 
unity  and  above,  eddy  formation  is  pronounced  and  a  scour 
hole  develops  completely  and  uniformly  about  the  pile  with 
a  maximum  depth  of  407o  of  the  pile  diameter. 

Palmer  noted  that  the  scour  hole  initiates  at  the 
sides  of  the  cylinder  and,  in  many  cases,  some  sediment 
deposition  was  found  to  occur  against  the  upstream  and  down- 
stream sides  of  the  cylinder.   As  the  scour  hole  continues 
to  develop,  these  deposits  are  removed  and  a  majority  of 
the  sediment  deposition  shifts  to  a  rim  outside  the  scour 
hole.   Once  fully  developed,  the  scour  hole  was  found  to 
be  a  radially  symmetric  frustrum  of  an  inverted  cone  hav- 
ing side  slope  at  the  angle  of  repose  of  the  sediment. 
Palmer's  data  indicated  that  the  depth  of  the  equilibrium 
scour  hole  to  the  pile  diameter  ranged  from  between  0.5 
and  0.75,  a  value  higher  than  that  previously  stated  by 
Ranc  e . 


For  oscillating  flow  past  a  cylindrical  obstruction, 
the  horseshoe  vortex  is  insignificant,  and,  therefore, 
scour  is  initiated  from  the  bottom  fluid  shear  stress 
caused  by  the  acceleration  of  primary  flow  past  the  ob- 
struction.  This  results  from  a  bottom  fluid  shear  stress 
gradient  that  is  positive  from  upstream  of  the  cylinder 
to  the  point  of  maximum  velocity,  and  negative  from  that 
point  to  the  downstream  limit  of  undisturbed  flow  (Figure 
10) .   If  the  maximum  bottom  fluid  shear  stresses  are  great- 
er than  the  sediment  entrainment  threshold,  erosion  will 
occur  from  upstream  of  the  cylinder  up  to  the  point  of 
maximum  velocity,  and  deposition  will  occur  downstream  of 
this  point. 

In  order  for  an  equilibrium  scour  hole  to  develop 
around  the  entire  circumference  of  the  cylinder,  an  addi- 
tional flow  condition  must  develop.   The  initial  scour  hole 
produces  bottom  slopes  close  to  the  pile  which  must  even- 
tually cause  flow  separation  in  the  thin  orbital  boundary 
layer.   Reverse  eddies  in  the  separated  flow  erode  the 
initial  side  slope  of  the  scour  hole  until  it  obtains  a 
value  of  the  angle  of  repose  of  the  sediment.   With  the 
steepening  of  the  sideslopes,  the  zone  of  separated  bottom 
boundary  layer  flow  spreads  out  around  the  pile  eventually 
causing  the  scour  hole  to  obtain  the  equilibrium  shape 
shown  in  Figure  11.   Therefore,  the  equilibrium  scour  hole 
for  oscillatory  flow  is  related  to  boundary  layer  separation 
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Fig,  10.  -  Sediment  transport  due  to  oscillatory  flow  (15) 
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Fig.    11.   -   Oscillatory   flow  scour   cavity    (15) 


and  not  dependent  on  the  nature  of  secondary  flow  as  was 
the  case  for  scour  due  to  steady  currents  (15) . 
Mechanics  of  Scour  Due  to  Waves  and  Currents 

In  this  study,  the  scour  that  developed  through  the 
combined  action  of  waves  and  currents  was  investigated. 
This  is  the  type  of  scour  that  occurs  in  almost  all  ocean 
environments.   Ninomiya  et  al.  (16)  performed  scale  model 
experiments  using  the  combined  action  of  waves  and  currents. 
They  found  the  shape  of  the  scour  hole  to  be  similar  in 
shape  to  that  developed  by  a  steady  current  without  waves, 
and  that  the  rate  of  scour  hole  development  was  consider- 
ably faster  in  the  case  of  combined  wave  and  current  action 
when  compared  to  steady  current. 

As  mentioned  earlier,  if  the  steady  flow  velocity  is 
stronger  than  the  oscillating  flow  velocities,  a  horseshoe 
vortex  will  form  about  the  base  of  the  pile.   The  intensity 
of  this  vortex  varies  according  to  differences  in  the 
square  of  the -combined  steady  and  oscillating  flow  velo- 
cities.  The  size  and  shape  of  the  equilibrium  scour  hole 
should  not  be  as  large  as  that  caused  by  a  steady  current 
whose  mean  velocity  is  equal  to  the  combined  velocity  of 
the  wave  orbital  and  mean  current  condition.   This  is  due 
to  the  fact  that  velocity  profile  of  the  combined  veloci- 
ties will  be  flatter  than  that  achieved  in  a  simple 
steady  boundary  layer  flow.   On  the  other  hand,  the  equi- 
librium scour  hole  depth  should  be  greater  than  that  caused 
by  a  steady  boundary  layer  having  no  wave  action. 
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If  a  strong  current  and  a  weak  wave  orbital  flow 
exists,  the  development  of  the  scour  hole  will  depend 
primarily  on  a  horseshoe  vortex  with  a  pulsing  intensity. 
The  scour  hole  development  will  be  faster  when  compared 
to  a  steady  current  alone  due  to  the  additional  bottom 
fluid  shear  stress  caused  by  the  wave  orbital  boundary 
layer  and  because  the  vortex  formed  from  the  steady  por- 
tion of  the  flow  can  more  easily  transport  the  wave- 
agitated  sediment. 

In  the  case  of  a  strong  orbital  flow  and  a  weak  cur- 
rent, the  scour  hole  shape  will  be  similar  to  that  caused 
by  a  steady  current.   Once  again,  the  scour  will  develop 
more  rapidly  than  it  would  under  the  influence  of  a  steady 
current  alone  due  to  the  additional  bottom  fluid  stresses 
caused  by  the  waves  orbital  motion.   In  this  case,  the 
equilibrium  scour  depth  can  be  less  than  that  associated 
with  a  steady  current  that  is  equal  to  the  sum  of  the  weak 
current  and  orbital  velocity. 

The  situation  becomes  difficult  to  evaluate  when  the 
magnitude  of  the  oscillatory  motion  and  the  steady  current 
are  nearly  equal.   The  scour  hole  size  and  shape  will  nor- 
mally depend  on  the  horseshoe  vortex  associated  with  a 
steady  current.   It  is  not  possible  to  generalize  on  the 
scale  and  intensity  of  the  pulsing  vortex  when  compared 
to  the  steady  current  situation. 
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The  size  and  depth  of  scour  produced  by  the  combination 
of  waves  and  currents  acting  on  a  pile  that  is  large  in 
comparison  to  the  oscillation  amplitude  will  be  drastical- 
ly smaller  than  that  predicted  for  a  small  pile,  except 
when  the  steady  current  is  greater  than  the  maximum  wave 
orbital  velocity.   This  is  caused  by  the  fact  that  the 
horseshoe  vortex  does  not  have  time  to  wrap  completely 
around  the  pile  before  the  flow  reverses  due  to  the  wave 
orbital  velocity  (15) . 
Previous  Studies  on  Scour 

Very  little  experimental  work  has  been  done  on  scour 
caused  by  the  combined  effect  of  oscillatory  wave  motion 
and  current;  however,  there  have  been  numerous  studies 
dealing  with  scour  in  open  channel  flow.   Because  the 
scouring  forces  are  similar  for  oscillatory  flow  as  for 
steady  state  open  channel  flow,  that  is,  hydrodynamic  in 
nature,  the  results  found  in  experiments  in  open  channel 
flow  can  be  applied  with  reservation  to  the  combined  oscil- 
latory and  steady  current  flow. 

Some  experimental  results  are  presented  by  Shen  et 
al .  (22)  as  shown  in  figures  12-14.   These  experiments 
employed  a  steady  current.   Figure  12  indicates  that  the 
depth  of  scour  is  a  function  of  the  pier  Reynolds  number 
(TR). 

TR  =  UD/v 
where 
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Fio.   i3#  _  Equilibrium  scour  depth  -  Pier  Reynolds 
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U  =  horizontal  free  stream  velocity 

D  =  pile  diameter 

v  =  kinematic  viscosity 
It  was  also  shown  that  the  scour  depth  is  also  dependent 
on  the  pile  diameter  (Figure  13)  and  sediment  size  (Figure 
14) ,  and  this  explains  why  the  Reynolds  number  alone  is 
insufficient  to  collapse  all  of  the  data  shown  in  Figure 
12  onto  a  single  line. 

Carstens  (3)  performed  studies  of  scour  associated 
with  various  types  of  obstacles.   He  concluded  that  the 
rate  of  scour  induced  by  the  obstacle  is  a  function  of 
the  sediment  number  (N  ) ,  sediment  geometry,  as  well  as 
the  ratio  of  the  scour  depth  to  obstacle  size.   It  should 
be  noted  that  the  "sediment  number"  is  a  densimetric  Froude 
numb  er . 


Ns  =  U/4   (Ss-l)gd 

where 

U  =  free  stream  velocity 
S  =  sediment  specific  gravity 

g  =  acceleration  of  gravity 

d  =  mean  diameter  of  sediment  particle 

Coleman  (5)  used  the  principle  that  in  any  problem 
of  flow  around  an  immersed  body,  where  the  dependent  vari- 
able is  a  pressure  drop,  the  most  general  formulation  of 
the  problem  is  in  terms  of  an  Euler  number  and  either  a 


JO 


Reynolds  number,  a  Froude  number  or  both,  with  the  dependent 
variable  appearing  in  the  Euler  number.   He  derived  a  rela- 


tion between  the  scour  Euler  number  (U/^l  2gd  )  and  the  pier 
Froude  number  (U/  V"gD)  . 

Zga~  -  °-6(^=>        (35) 

&  s        VgD 
From  the  above  relation  and  using  a  previously  derived 
depth  of  scour  relation  using  the  pier  Froude  number,  an 
expression  was  obtained  to  predict  the  depth  of  scour  for 
steady  flow. 

d  =  1.49  D9/1°  (U2/2g)1/10 (36) 

where 

d  =  depth  of  scour  (ft) 

D  =  pile  diameter 

U  =  free  stream  velocity 

g  =  acceleration  of  gravity 
Wells  and  Sorensen  (27)  investigated  scour  around  a 
single  circular  pile  due  to  oscillatory  wave  motion.   Fig- 
ure 15  shows  that  a  relationship  exists  between  the  rela- 
tive scour  and  the  relative  water  depth  for  various  sedi- 
ments.  In  addition,  their  results  indicated  that  the  ra- 
tio of  ultimate  scour  to  wave  height  appeared  to  be  di- 
rectly related  to  the  sediment  number,  NgI  (Figure  16)  and 
to  the  pile  Reynolds  number,  NRp ,  (Figure  17).   In  Figures 
15  though  17,  the  notation  is  defined  as  follows: 
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Su  -  ultimate  significant  scour  depth 
d  =  mean  grain  diameter 

Sand  no .  1:   d  =  . 62  mm 
Sand  no .  2:   d  =  .325  mm 
Sand  no  .  3 :   d  =  .30  mm 

h  =  still  water  depth 
g  =  acceleration  of  gravity 
T  =  wave  period 
H  =  wave  height 

U  =  maximum  free  stream  velocity 
S   =  specific  gravity  of  sand 
D  =  pile  diameter 
v  =  kinematic  viscosity  of  water 
More  recent  scour  studies  include  the  work  of 
Raudkivi  and  Ettema  (20)  ,  Jain  and  Fischer  (11) ,  Armbrust 
(2) ,  and  Wang  (26)  .   The  work  of  Armbrust  and  Wang  dealt 
with  the  scour  associated  with  the  combined  effect  of  waves 
and  currents  while  the  remaining  studies  dealt  only  with 
a  steady  current. 

Raudkivi  and  Ettema  (20)  performed  laboratory  experi- 
ments on  clear  water  scour  of  cohesionless  bed  sediment 
at  cylindrical  piles.   They  found  that: 

1.   As  the  geometric  standard  deviation  of  the  par- 
ticle size  distribution  increases  then  the  equi- 
librium scour  depth  decreases. 


2.  The  equilibrium  depth  of  scour  decreases  with 
the  decreasing  relative  size  of  pier  diameter 
to  mean  bed  particle  size  for  values  less  than 
about  20  to  25. 

3.  The  maximum  equilibrium  depth  of  scour  is  a 
function  of  the  relative  depth  h/D  and  the  rela- 
tive pier  size  D/D,-n. 

4.  The  equilibrium  depth  of  scour  decreases  at  a 
greater  rate  with  decreasing  flow  depth  for 
smaller  values  of  the  relative  flow  depth  h/D. 

Jain  and  Fischer  (11)  performed  scour  experiments 
for  steady  currents  at  high  velocities,  i.e.,  for  flow 
regimes  well  beyond  the  onset  of  general  sediment  trans- 
port.  From  these  studies,  they  concluded  that  after  de- 
creasing slightly  upon  the  start  of  sediment  transport, 
scour  depth  does  increase  again  with  increasing  velocities 
This  conflicts  with  the  consensus  of  previous  experiments 
which  concluded  that  once  the  conditions  of  general  sedi- 
ment motion  have  been  established  in  a  stream,  there  will 
be  no  further  increase  of  scour  depth  with  velocity,  be- 
cause the  rate  at  which  sediment  enters  the  scour  region 
is  equal  to  the  rate  of  sediment  removal  by  the  scouring 
process . 

Armbrust  (2)  found  that  in  general,  for  a  given  cur- 
rent and  relative  depth  (h/L) ,  scour  increases  as  the 
wave  steepness  (H/L)  increases  (Figure  18) .   Figure  18 
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also  shows  that  the  ultimate  scour  depth  increases  as  the 

relative  depth  decreases.   Additionally,  in  general,  for 

a  fixed  relative  depth,  Figure  19  reveals  an  increase  of 

scour  depth  with  an  increase  of  relative  wave  height  (H/h) . 

In  Figures  18  and  19,  the  notation  is  defined  as  follows: 

S   =  ultimate  scour  depth 
u  r 

D  =  pile  diameter 
Dcq  =  mean  sediment  diameter 
H  =  wave  height 
L  =  wave  length 
U-,  ^  =  free  stream  velocities 
h  -  still  water  depth 
Wang  (26)  defined  ultimate  scour  (S  )  to  include 
"general  scour"  and  net  "local  scour"  and,  therefore,  ul- 
timate scour  is  equal  to  the  total  scour  that  occurs.   Wang 
stated  that: 

Su  =  f(Vfw,h,p,H,L,T,d50,Ps,D,v,t) 

where 

Vr   =  combined  wave  and  current  velocity 
fw 

h  =  water  depth 

p  =  water  density 

H  =  wave  height  with  current 

L  =  wave  length  with  current 

T  =  wave  period  with  current 

D  =  pile  diameter 

dcn  =  mean  sediment  diameter 
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p   =  density  of  sediment 

s         J 

v  =  kinematic  viscosity  of  water 
t  =  time  required  to  produce  scour 
Using  dimensional  analysis,  Wang  produced  the  relation 


S 

-£  =  f (NF  •  H/L  •  Ur  •  Ng  •  NRp)  =  f (8) 


where 


=  NF  •  H/L  •  Ur  •  Ns  •  Nrp 

V3H2L[V  +  (1/T  -  V/L)  ^-]2  D 

Ps  -  '   2  a        ^ 
n §  vdsn  h 
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Figure  20  is  a  plot  of  the  relationship  of  relative  scour 


depth,  S  /h,  and  the  parameter 
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PART  III 
THEORETICAL  ASPECTS 
Wave-Current  Interaction 

Hoften  and  Karaki  (23)  performed  a  laboratory  study- 
to  investigate  the  interaction  of  waves  and  currents  as 
mechanically  generated,  monochromatic  waves  propagated 
on  a  turbulent  open  channel  flow. 

They  found  that  waves  propagating  on  an  open  channel 
flow  are  shown  to  be  modified  by  a  combination  of  the  wave 
interaction  with  the  shear  gradient  and  the  existing  tur- 
bulence.  The  effect  was  to  remove  energy  from  the  waves 
causing  an  increase  in  wave  energy  dissipation. 

The  flow  itself  is  also  changed  by  the  waves.   The 
induced  orbital  wave  velocities  distort  the  mean  velocity 
profile  increasing  the  gradient  close  to  the  sediment  bed. 
Turbulence  energy  is  formed  near  the  surface  by  the  in- 
teraction of  the  wave -induced  Reynolds  stresses  and  the 
fluctuating  velocity  gradients.   The  energy  diffuses  down- 
ward where  it  is  dissipated  on  the  sediment  bed  in  the 
form  of  an  apparent  increase  in  boundary  shear  stress,  an 
important  factor  in  the  scouring  process. 

In  order  to  evaluate  the  combined  velocity  of  a  wave 
propagating  in  the  same  direction  as  a  uniform  current, 
the  Airy  wave  theory  will  be  used.   By  doing  so,  if  a 
uniform  current  is  assumed,  the  velocity  potential  of  the 
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uniform  current  (<j>,)  and  the  velocity  potential  due  to 
wave  (tfO  can  be  added  directly  to  obtain  the  velocity 
potential  of  the  combined  wave  current  condition  (<$>)    be- 
cause both  velocity  potentials  are  linear  solutions  of  the 
Laplace  equation.   The  parameters  of  wave  period  (T) ,  wave 
height  (H)  and  wave  length  (L)  are  shown  in  Figure  (21) , 
and  are  properties  of  a  wave  propagating  in  still  water. 
Wave  velocities  can  be  calculated  using  these  parameters 
with  the  water  depth  for  a  no  current  condition.   There- 
fore, the  velocity  relation  for  combined  wave-current  in- 
teraction becomes, 

♦  -  *i  +  *2 (12) 

where 

<j>,  =  U  x (13) 

Yl    c 

and  the  derivation  of  <t> «  is  outlined  in  the  following  pages 
This  procedure  closely  follows  Newman  (14) . 

The  fluid  velocity  under  wave  action  is  expressed  by 
the  gradient  of  the  velocity  potential  <j>  ,  and  the  effects 
of  the  free  surface  must  be  expressed  in  terms  of  appropri- 
ate boundary  condition  on  this  surface. 

The  exact  kinematic  boundary  condition  can  be  derived 
most  readily  by  requiring  that  the  substantial  derivative 
of  the  quantity  z  -  n  vanish  on  the  free  surface.   The 
result  of  this  condition  is  that  on  z  =  n, 

D      /             n_34>        f^L       i^ilL       jty.£jL  (14) 

0   "   Dt    (z   "    n)         "3z    "    9t    "    8x    3x         ay  ay 
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For  the  linearized  kinematic  boundary  condition,  the  second 
order  terms  can  be  dropped,  which  gives  the  relation: 

df]  3(f) 

3t  ~  Si" (15) 

This  approximate  boundary  condition  states  that  the  verti- 
cal velocities  at  the  free  surface  and  the  fluid  particles 
are  equal,  disregarding  small  departures  of  the  surface 
from  the  horizontal  orientation. 

The  dynamic  boundary  condition  is  obtained  from 
Bernoulli's  equation.   The  exact  condition  to  be  satis- 
fied on  the  free  surface  is: 

-  i  (P  -  Pa)  =  |1  +  1/2  v  *  .  v  *  +  gz  =  0 (16) 

Substituting  the  free  surface  elevation  n  for  Z  and  solv- 
ing for  n  the  expression  becomes, 

n  =  -1/g  (||-+  1/2  7  ♦  .  v  *) (17) 

By  neglecting  the  last  term,  which  is  of  the  second  order, 
the  linearized  equation  for  the  free  surface  elevation  is 

"  =  -l/S  H  • (18) 

On  the  surface,  Z  =  0,  the  linearized  dynamic  boun- 
dary condition  (18)  can  be  differentiated  with  respect  to 
time  and  added  to  the  linearized  kinematic  boundary  con- 
dition (15) ,  to  form  a  single  boundary  condition  for  the 
velocity  potential: 


Jl 


324>  d<t>  n 

— T+  §    ai"  =    0   on   z   =    0 (19) 

dt 

For  plane  progressive  waves  with  finite-depth  effects, 
the  free  surface  elevation  will  be  of  the  general  form, 

n(x,t)  =  a  cos  (kx  -  oat  +  e)   (20) 

in  which  the  positive  x  -  axis  is  in  the  direction  of  wave 
propagation,  a  is  the  wave  amplitude,  and  e  is  some  arbi- 
trary phase  angle  that  can  set  equal  to  zero  by  the  suit- 
able choice  of  the  origin  -x  =  0.   The  parameter  k  is  the 
wave  number,  in  which, 

k  =  2tt/l 

where 

L  =  wavelength 

The  solution  to  this  problem  is  written  in  terms  of 
a  two  dimensional  velocity  potential  <f>(x,z,t),  which  must 
satisfy  the  free-surface  condition  (19),  Laplace's  equa- 
tion, the  bottom  boundary  condition  3^/8 z  =  0  at  z  =  -d 
(for  finite  depth) ,  and  must  produce  the  wave  elevation 
(20)  for  equation  (18) .   The  potential  must  be  sinusoidal 
due  to  (20)  and  the  solution  should  be  in  the  form: 

4>(x,z,t)  =  Re[Z(z)e'ikx  +  iwt]   (21) 

From  Laplace's  equation,  Z  must  satisfy  the  ordinary  dif- 
ferential equation: 

1-f  -  k2Z  =  0 (22) 

dzz 
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The  most  general  solution  to  equation  (22)  in  terms 
of  exponential  functions  is 

Z  =  Cekz  +  De"kz (23) 

where 

C  =  constant 

D  =  constant 

First  assume  the  fluid  depth  is  infinite  or 
-°°  <  z  <  0.   To  avoid  an  unbounded  motion  deep  below  the 
free  surface,  the  constant  D  in  equation  (23)  must  vanish 
leaving 

Z  =  Cekz • (24) 

If  equation  (24)  is  substituted  into  equation  (21) ,  the 
velocity  potential  becomes, 

*  =  Re[Cek2"ikx+iwt]  (25) 

By  setting  z  =  0,  using  equation  (18)  to  calculate  the 
linearized  wave  elevation,  and  comparing  the  result  with 
the  required  plane  wave  in  equation  (20)  ,  the  constant 
C  has  the  value  of 

C  =  iga/oa (26) 

For  a  fluid  having  a  constant  depth  d,  with  the  bot- 
tom being  an  impermeable  plane,  as  stated  earlier,  the 
bottom  boundary  condition  becomes, 


|£  =  0   at  z  =  -d (27) 

By  looking  at  the  general  solution  (equation  (23)),  it  is 
found  that  both  exponential  functions  must  be  retained, 
with  constants  C  and  D  chosen  such  that  equation  (27)  will 
be  satisfied. 

For  simplification,  the  two  independent  solutions 

+kz 

e  —   are  replaced  by  linear  combinations  of  the  hyperbolic 

sine  and  cosine  function.   For  a  fluid  having  finite  depth, 
the  solution  becomes, 

m        cosh  k  (z  +  d) '   (28) 

v  '       cosh  kd 

The  constant  C  is  selected  such  that  equation  (26)  remains 
valid.   When  substituting  (28)  into  (21)  and  using  equa- 
tion (26)  the  velocity  potential  for  a  finite  depth  fluid 
becomes , 

m   ga  cosh  k(z  +  d)  s.n  (kx  _  ut)   (29) 

^2    co     cosh  kd 

Substituting  equation  (29)  and  (13)  into  equation 
(12),  the  combined  wave-current  velocity  potential  becomes: 

^  =  Ux  +  ^  cosh  k(z  +  d)  sin  (kx  _  ut) (30) 

*    c     u     cosh  kd 

Differentiating  equation  (30),  the  result  becomes: 

n  =  _3±  =  n  +  ^  cosh  k(z  +  d)   Cos  (kx  -  a)t)  .  .  .  .  (31) 
u    3x    cm       cosh  kd 

where 


JM- 


U  =  combined  wave-current  velocity 

It  can,  therefore,  be  concluded  that  the  combined 
wave-current  velocity  (U)  can  be  obtained  by  simply  add- 
ing the  value  for  the  uniform  current  velocity  to  the 
wave  velocity  in  still  water  found  using  the  parameters 
stated  earlier. 

It  is  appropriate  at  this  time  to  mention  several 
points.   Firstly,  this  relationship  is  invalid  for  non- 
linear boundary  conditions.   Also,  since  the  wave  pro- 
pagates into  the  flow,  there  will  be  modifications  to 
both  the  waves  and  currents  by  interaction,  therefore, 
the  simple  addition  of  current  velocity  and  wave  velocity 
is  only  approximate.   Finally,  linear  wave  theory  is  only 
valid  for  small  amplitude  waves,  and  many  cases  of  wave- 
current  interaction  exist  beyond  the  accepted  range  of 
validity  for  this  theory. 

Wang  (26)  calculated  the  combined  wave-current  velo- 
city using  the  procedure  outlined  above  and  compared  the 
calculated  values  with  measured  values .   He  derived  a 
combined  wave-current  velocity  relation  given  by: 

U  =  Uc  +  ttHU/T  -  Uc/L)  slnh121  d   cosh[^  (z  +  d)]   (32) 

L 
where 

U  =  combined  wave-current  velocity 

U  =  uniform  current  velocity 
c 

By  assuming  that  the  distribution  of  current  is  uni- 
form with  depth,  the  depth-averaged  velocity  was  obtained 


from  the  integral  of  equation  (32)  .   When  comparing  the 
calculated  and  measured  combined  wave-current  velocities, 
a  correlation  coefficient  of  0.939  was  obtained,  thus  in- 
dicating this  method  for  calculating  currents  is  reasonable 

Dimensional  Analysis 

There  are  many  variables  encountered  when  studying 
scour.   The  following  list  is  believed  to  be  the  important 
variables  influencing  scour  depth  at  a  cylindrical  pile 
due  to  current  and  random  wave  motion. 


Description 
Water  depth 


Symbol 

Units 

Dimensions 

h 

ft 

L 

Hl/3 

ft 

L 

V 

ft2/s 

L2/T 

D 

ft 

L 

D50 

ft 

L 

P 

3 

slugs/ft 

M/L3 

(Ps  -  P) 

3 

slugs/ft 

M/L3 

g 

ft/sec 

L/T2 

U 

ft/sec 

L/T 

Sj 

ft 

L 

Significant  wave 
height 

Kinematic  viscosity 
of  water 

Pile  diameter 

Mean  diameter 
of  sediment 

Density  of  water 

Sediment/ fluid 
density  difference 

Acceleration  due 
to  gravity 

Combined  wave- 
current  velocity 

Ultimate  scour 
depth 

Since  random  waves  were  used  in  this  study,  and  they 
are  defined  by  the  significant  wave  height  (H-^3) ,  wave 


DO 


period  and  wave  length  have  been  omitted  from  the  above 
list  of  variables. 

The  dimensional  analysis  was  performed  using  a  me- 
thod outlined  by  Venezian  (25)  during  the  1984  Texas  A&M 

University  dredging  short  course.   S,  was  chosen  as  the 

2  3 

unit  of  length,  S,  /v  as  the  unit  of  time,  and  p/S,   as 

the  unit  of  mass.   The  remaining  seven  variables  expressed 
in  these  units  give  the  following  dimensional  numbers: 

Sd/h,  Sd/H1/3,  Sd/D,  Sd/D5Q,  p/(ps  -  p),  v2/Sd3g,  v/Sdu 

By  rearranging  and  combining  terms : 

(Sd/D)  (v/Sdu)    v    LNRP 
[(Sd/h)(l/(v/Sdu)2)(v2/Sd3g)]%  =  u/^fgh   =  NF 

[(Sd/D)(l/(v/Sdu)2)(v2/Sd3g)]%  =  u/Yp  =  NFp 

[(%)(v2/Sd3g)(l/(v/Sdu)2)]%    =  u/V2S^i  =  NES 

[(l/(v/Sdu)2)(v2/Sd3g)(Sd/D50)(P/(ps  -  p))]% 


=  U/VgD50((Ps  "  p)/p)  =  Ns 
Therefore,  the  functional  relationship  for  scour  is: 

TT=  f(Sd/Hl/3>  Sd/D>  Sd/D50>  NRP'  V  NFP>  NES'  V 


3/ 


PART  IV 
EXPERIMENTAL  EQUIPMENT 

The  scour  experiments  were  conducted  in  a  120-foot 
long,  3-foot  deep  and  2-foot  wide  two  dimensional,  glass- 
walled  wave  flume.   The  test  area  was  located  approximate- 
ly 40  feet  from  the  wave  generator.   The  entire  test  area 
measured  40  feet  in  length  and  was  divided  into  several 
sections.   An  8-foot  long  5-inch  deep  test  sediment  bed 
was  bordered  on  each  end  by  an  8-foot  long  10- inch  high 
false  bottom.   At  the  leading  edge  of  the  test  area,  a 
16-foot  long  ramp  was  used  to  gradually  bring  the  wave 
up  to  the  new  depth  at  the  top  of  the  false  bottom.   A 
schematic  of  the  test  area  is  shown  in  Figure  (22) .   A 
1.5-inch  aluminum  pile  was  positioned  along  the  center- 
line  of  the  sediment  bed. 

Currents  were  produced  using  a  centrifugal  pumping 
system.   The  flow  rate  was  monitored  using  a  Fischer  and 
Porter  electromagnetic  flow  meter.   The  pump  discharge 
entered  the  wave  flume  just  forward  of  the  wave  generator, 
and  was  diffused  using  a  louvered  grate.   A  Thermo -Systems 
model  1050  hot  film  anemometer  was  used  to  measure  the 
vertical  velocity  distribution  of  the  combined  waves 
and  currents . 

The  random  wave  generator  used  was  a  Seasim  modular 
wave  making  system.   This  system  consisted  of  a  servo 
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control  system  amplifier,  a  Seasim  programmable  spectrum 
random  signal  generator,  a  Seasim  autocompensating  wave 
height  gauge,  and  a  paddle  type  sea  wave  simulation  roll- 
ing seal  modular  wave  maker.   A  capacitance-type  wave  gauge 
was  used  with  a  Hewlett-Packard  model  17403A  carrier  pre- 
amplifier to  measure  the  wave  profiles.   To  determine  the 
significant  wave  height  from  the  wave  profile,  a  wave 
digitizer  was  utilized. 

A  wave  absorber  was  positioned  in  the  wave  flume  at 
the  opposite  end  of  the  wave  generator  to  reduce  any  re- 
flections which  might  occur. 

Scour  depth  readings  were  made  using  a  point  gauge 
after  each  test. 


ou 


PART  V 
EXPERIMENTAL  PROCEDURES 

A  total  of  20  experiments  were  performed  during  this 
investigation.   Four  runs  utilized  currents  only,  and  the 
remaining  16  runs  used  a  combination  of  waves  and  cur- 
rents.  All  of  the  test  runs  used  a  1.5-inch  diameter  pile 
placed  in  a  bed  of  glass  microbeads .   The  glass  microbeads 
had  a  mean  diameter  (D,«)  of  0.1  mm  and  a  specific  gravity 
of  2.45.   For  both  the  current  only  tests  and  the  combined 
wave-current  tests  two  water  depths  and  two  different  cur- 
rents at  each  depth  were  used.   Also  used  were  four  dif- 
ferent random  wave  spectra  for  each  depth/current  combina- 
tion.  The  data  on  the  four  random  wave  spectra  are  con- 
tained in  Appendix  A,  and  Table  1  contains  a  description 
of  all  tests  conducted. 

At  the  commencement  of  each  test,  the  sediment  bed 
was  leveled  to  a  5-inch  uniform  thickness.   Water  was 
added  slowly  to  the  tank  so  as  not  to  disturb  the  sedi- 
ment.  After  the  required  water  depth  was  attained,  the 
proper  flow  rate  through  the  tank  was  established.   The 
water  level  in  the  tank  was  maintained  by  regulating  a  10 
inch  discharge  valve  located  at  the  rear  of  the  tank. 
Once  the  proper  flow  rate  was  attained  and  the  required 
water  depth  held  constant,  the  sediment  bed  was  once  again 
leveled.   The  signal  generator,  which  had  been  previously 
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TABLE  1  -  SUMMARY  OF  TESTS 


Test 

Current 

1 

ui 

2 

U2 

3 

U3 

4 

U4 

5 

Ul 

6 

U2 

7 

U2 

8 

U2 

9 

U2 

10 

Ul 

11 

Ul 

12 

Ul 

13 

U4 

14 

U4 

15 

U4 

16 

U4 

17 

U3 

18 

U3 

19 

U3 

20 

U3 

Ul  - 

.405  fps 

U2  =  .380  fps 
U3  =  .650  fps 
U4  =  .485  fps 

D  =  Darbyshire 

P  =  Pierson  -  Moskowitz 
Jl  =  Jonswap  (T  =  8  sec.) 
J2  =  Jonswap  (T  =  10  sec.) 


h1  =  1.167  ft 
h2  =   .917  ft 


Wave  Spectrum         Depth 

hl 
h2 

h2 

hl 
D  hx 

D  h2 

P  h2 

Jl  h£ 

J2  h2 

P  h. 


Ll 
J2  h. 


Jl  hx 


ll 

D  h1 

P  h, 

Jl  \ 

J2  h1 

D  h2 

P  h2 

Jl  h2 

J2  h2 


o± 


programmed,  was  started  generating  the  desired  wave  spec- 
trum.  The  maximum  scour  depth  was  recorded  at  various  in- 
tervals throughout  the  run,  the  time  increment  increasing 
as  the  test  continued.   Each  test  was  conducted  until  an 
equilibrium  in  the  measured  maximum  scour  depth  was  at- 
tained.  Following  the  completion  of  the  test,  the  wave 
generator,  and  the  flow  through  the  tank  were  stopped  and 
the  water  was  slowly  drained  from  the  tank  in  order  to 
preserve  the  scour  pattern  formed.   Once  the  tank  was  com- 
pletely drained,  photographs  were  taken  of  the  scour  pat- 
terns .   In  addition,  sediment  elevations,  taken  in  the  form 
of  a  rectangular  grid,  were  measured  using  a  point  gauge 
in  order  to  determine  the  equilibrium  scour  depth.   This 
depth  was  used  in  the  data  calculations  presented  later 
in  this  paper. 


PART  VI 
PRESENTATION  AND  DISCUSSION  OF  RESULTS 

The  experimental  data  obtained  from  this  study  are 
presented  below  in  Table  2: 

TABLE  2  -  EXPERIMENTAL  DATA 

Run      Velocity  (fps)      Scour  Depth  (ft)     Hl/3  (^ft') 

1  U  =  .405  .101 
c 

2  U  =  .380  .081 
c 

3  Uc  =  .650  .167 

4  Uc  =  .485  .14 

5  U  =  .442  .126  .2278 

6  U  =  .394  .115  .1656 

7  U  =  .393  .13  .1266 

8  U  =  .401  .127  .1715 

9  U  =  .386  .09  .1930 


Velocity    (fps) 

U 
c 

U 
c 

= 

.405 
.380 

U 
c 

= 

.650 

U 
c 

= 

.485 

U 

= 

.442 

u 

= 

.394 

u 

= 

.393 

u 

= 

.401 

u 

=r 

.386 

u 

= 

.437 

u 

= 

.459 

u 

= 

.419 

u 

= 

.498 

u 

= 

.491 

u 

= 

.504 

u 

= 

.508 

u 

= 

.665 

u 

= 

.659 

u 

= 

.660 

u 

= 

.664 

10  U  =  .437  .122            .1817 

11  U  =  .459  .152            .1791 

12  U  =  .419  .113            .2100 

13  U  =  .498  .143            .1724 

14  U  =  .491  .149            .1527 

15  U  =  .504  .174            .1631 

16  U  =  .508  .155            .1550 

17  U  =  .665  .18             .1342 

18  U  =  .659  .174            .0967 

19  U  =  .660  .176            .1605 

20  U  =  .664  .168            .1093 

The  results  of  the  laboratory  experiments  are  present- 
ed in  graphical  form  in  this  section.   These  graphs  depict 
the  relationship  between  the  dimensionless  parameters 


developed  earlier  for  scour  as  well  as  several  dimensional 
plots.   On  various  dimensionless  plots,  data  from  experi- 
ments conducted  by  Wang  (22)  and  Armbrust  (2)  were  used 
in  combination  with  data  obtained  from  the  random  wave- 
current  interaction  experiment.   On  several  plots,  curves 
were  constructed  to  indicate  a  trend  only,  and  these  plots 
were  labeled  accordingly.   Regression  analysis  was  per- 
formed on  the  remaining  plots  in  order  to  construct  the 
appropriate  curves.   Logarithmic  regression  analysis  was 

used  on  the  semi-logarithmic  plots,  and  power  regression 

2 
analysis  was  used  on  the  full  logarithmic  plots.   The  r  " 

term  adjacent  to  each  curve  represents  the  correlation 

coefficient  for  each  curve. 

Figures  (23)  to  (26)  show  the  functional  relationship 
that  exists  between  the  relative  scour  depth,  Sj/h,  and 
Sj/D  and  the  sediment  number  for  three  distinct  pile  dia- 
meters and  three  different  sediment  sizes.   It  shows  that 
scour  depth  increases  with  increasing  wave-current  velo- 
cities.  Scour  depth  is  not  dependent  on  the  sediment  num- 
ber alone,  but  is  also  dependent  on  pile  diameter  (Figures 
23  and  25)  and  sediment  size  (Figures  24  and  26) . 

Since  the  basic  mechanism  for  scour  hole  development 
is  the  vortex  system  resulting  from  flow  separation,  and 
the  separation  point  is  dependent  on  the  Reynolds  number, 
the  strength  of  the  vortex  system  is  a  function  of  the 
pile  Reynolds  number.   Because  of  the  relationship  between 
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the  vortex  system  and  the  pile  Reynolds  number,  scour  depth 
should,  in  turn,  be  functionally  related  to  the  pile  Rey- 
nolds number.   This  functional  relationship,  along  with 
Shen's  curve  for  scour  due  to  a  steady  current,  is  de- 
picted in  Figure  (27) . 

Figures  (28)  through  (31)  indicate  the  relationship 
between  relative  scour  depth,  S  ,/D  and  S,/h,  and  the  pile 
Reynolds  number.   It  shows  that  scour  depth  increases  with 
increasing  wave-current  velocities  or  pile  diameter  as  in 
the  case  of  the  sediment  number;  it  can  be  seen  from  the 
graphs  that  scour  depth  is  not  dependent  on  the  pile 
Reynolds  number  alone,  but  is  also  dependent  on  the  pile 
diameter  (Figures  28  and  30)  and  the  sediment  size  (Figures 
29  and  31) .   This  agrees  with  the  experimental  work  of 
Shen  et  al .  (22);  however,  the  dependence  of  the  pile 
Reynolds  number  on  sediment  size  differs  from  the  conclu- 
sions reached  by  Palmer  (18)  in  which  he  stated  that  the 
scour  depth  was  independent  of  grain  sizes  from  0.12  mm 
to  0.63  mm. 

By  comparing  Figures  (23)  to  (26)  and  Figures  (28)  to 
(31)  with  figures  found  in  Lee  (13)  which  plotted  the  same 
variables  for  multiple  piles,  it  is  noted  that  the  single 
pile  condition  requires  a  higher  velocity  to  produce  the 
same  scour  depth  as  the  multiple  pile  case,  assuming  all 
other  variables  (i.e.,  D5Q,h)  are  equal. 
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It  was  found  that  if  the  pile  Reynolds  number  and 
the  sediment  number  were  multiplied  together  and  plotted 
against  the  relative  scour  depth,  S,/h,  all  the  points 
collapsed  onto  a  single  curve  regardless  of  variance  in 
pile  diameter  or  mean  sediment  size.   This  result  is 
shown  in  Figures  (32) ,  (33)  and  (34) .   The  equation  for 
each  curve  is  also  indicated  on  each  figure.   Due  to  the 
tendency  of  the  data  for  random  waves  and  currents  and 
monochromatic  waves  and  currents  to  collapse  onto  a  sin- 
gle curve,  it  appears  as  if  the  irregularity  of  the  waves 
is  insignificant  in  terms  of  scour  depth  prediction.   By 
using  a  dimensionless  plot,  as  depicted  in  Figures  (32)- 
(34) ,  it  appears  as  though  predictions  of  scour  depth 
could  be  made  for  situations  in  which  the  pile  diameter 
or  mean  sediment  size  varied  without  utilizing  scaling 
factors.   Before  any  conclusions  can  be  made,  field  data 
are  required  to  evaluate  the  need,  if  any,  for  scale  fac- 
tors . 

Figures  (35)  through  (37)  show  the  relationship  found 
for  random  wave-current  interaction  between  the  relative 
scour  depth  S,/H-,,^  and  the  sediment  number,  N  ,  pile 
Reynolds  number,  N„p ,  and  (N  ) (N^p) ,  respectively.   The 
curves  indicate  that  there  is  a  rapid  initial  increase  of 
the  relative  scour  from  the  point  of  incipient  motion  up 
to  a  maximum  relative  scour  depth.   A  further  increase  of 
N  ,  NDT3  or  (N  )  (NDT3)  will  cause  a  rapid  decrease  in  the 

S     K.T  S     K-i: 
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Fig.  32.  -  Relative  scour  depth,  S,/h,  vs.  sediment 

number  times  pile  Reynolds  number  for  a  constant 
pile  diameter  and  various  sediment  sizes 
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relative  scour  depth,  this  decrease  being  independent  of 

Ns'  Nrp  or  ^Ns^NRp)  '   These  results  are  similar  to  those 
found  by  Wells  and  Sorenson  (27)  and  presented  earlier  in 
this  report  in  Figures  (16)  and  (17) .   The  relative  scour 
depth  is  smaller  in  Wells  and  Sorenson' s  case,  but  this 
can  be  explained  by  the  fact  that  they  used  wave  action 
only  and,  therefore,  their  scour  velocities  were  smaller 
in  magnitude.   The  pile  Reynolds  number  for  the  rapid  de- 
crease  of  relative  scour  depth  is  approximately  7.5  x  10 
in  both  cases,  this  rapid  decrease  in  S,/H-,yo  being  indepen- 
dent of  the  pile  Reynolds  number  and  the  sediment  size. 

A  dimensional  plot  was  constructed  in  Figure  (38) 
showing  the  relationship  of  ultimate  scour  depth,  Sj,  ver- 
sus the  combined  wave-current  velocity,  U.   It  was  found 
that  in  all  cases  the  equilibrium  scour  hole  was  deeper 
than  that  associated  with  a  steady  boundary  layer  current 
with  the  same  mean  velocity.   These  results  agree  with 
predictions  made  by  Niedoroda  and  Dalton  (15) .   In  gener- 
al ,  it  can  be  stated  that  the  scour  depth  caused  by  wave- 
current  interaction  is  approximately  107o  greater  than  the 
scour  depth  due  to  a  steady  current  alone. 

In  Figure  (39),  scour  depth  was  compared  to  H-,  ,~. 
Although  the  data  are  limited,  it  appears  as  though  for 
current  velocities  greater  than  0.45  feet  per  second,  there 
is  a  tendency  for  scour  depth  to  increase  with  increasing 
wave  height,  for  velocities  less  than  0.45  feet  per  second, 
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scour  depth  decreases  with  increasing  wave  height  and  for 
velocities  approximately  equal  to  0.45  feet  per  second, 
there  is  no  definite  trend  at  all.   This  result  is  most 
likely  due  to  the  fact  that  waves  and  currents  interact 
differently  as  the  magnitude  of  the  current  changes.   Al- 
though a  discussion  was  presented  earlier  in  this  paper 
on  the  mechanics  of  scour  due  to  waves  and  currents,  there 
is  still  a  great  deal  of  uncertainty  in  this  area,  and 
further  research  is  indicated. 

Scour  Pattern  Observations 

The  rate  of  scour  hole  development  by  wave-current 
interaction  was  faster  than  that  due  to  a  steady  current 
alone.   A  typical  example  of  scour  development  is  shown 
in  Figure  (40) .   As  seen  in  this  figure,  scour  development 
is  very  rapid  initially,  the  scour  rate  then  decreases 
until  the  equilibrium  scour  depth  is  reached. 

In  all  cases,  the  scour  hole  due  to  wave-current  in- 
teraction was  similar  in  shape  to  the  scour  hole  formed  by 
a  steady  current.   This  shape  resembles  an  inverted  cone. 
A  typical  example  of  a  scour  hole  formed  by  a  steady  cur- 
rent is  shown  in  Figure  (41) .   Figures  (42)  to  (45)  depict 
the  scour  pattern  formed  by  the  same  current  used  in  Run 
4  (Figure  41) ,  combined  with  four  different  random  wave 
spectra.   Photographs  of  the  scour  patterns  formed  by 
wave-current  interactions  for  the  remaining  12  runs  can 
be  found  in  Appendix  B.   As  stated  earlier,  the  equilibrium 
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RUM  4 


Fig,  Zf1.  -  Typical  scour  cavity  due  to  a 
steady  current 


h  =  1 .167  ft 
U  =  .485  ft/sec 
D=  .125  ft 
Sd=  .14  ft 
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Fig.    Zf2.   -   Scour  pattern  due   to  random  waves  and 
current  -  Run  13. 
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ig.   k3.   -  Scour  pattern  due   to  random  waves  and  current 
Run   1/f. 
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Fig.  Zf5.  -  Scour  pattern  due  to  random  waves 
and  current  -  Run  16. 
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scour  hole  associated  with  the  wave-current  interaction 
was  deeper  when  compared  to  the  equilibrium  scour  hole 
formed  by  a  steady  current  alone.   The  size  and  shape  of 
the  equilibrium  scour  hole  associated  with  wave-current 
interaction  was  at  times  greater  than,  approximately  equal 
to,  or  less  than  the  size  and  shape  of  the  scour  hole  due 
to  a  steady  current  alone.   These  results  differ  from  the 
prediction  made  by  Niedoroda  and  Dalton  (15) .   They  stated 
that  the  scour  hole  size  and  shape  for  combined  waves  and 
currents  should  not  be  as  great  as  that  developed  by  a 
steady  current  alone. 

Scour  development  is  illustrated  by  Figures  (46)  to 
(48) .   Figure  (46)  shows  a  typical  pattern  of  scour  develop' 
ment  in  its  initial  stages,  approximately  the  first  2  to 
5  minutes  of  the  test.   Figure  (47)  and  (48)  are  the  typi- 
cal scour  patterns  that  occur  due  to  wave-current  interac- 
tion 30  minutes  and  2  to  3  hours  after  test  commencement, 
respectively . 

The  general  scour  pattern  was  one  of  localized  scour 
around  the  pile.   A  vortex  formed  around  the  pile  due  to 
disturbed  flow  and  suspended  the  sediment.   The  sediment 
was  carried  farther  downstream  under  the  combined  effects 
of  waves  and  currents  than  it  was  with  a  steady  current 
alone.   Initially,  the  sediment  suspension  appeared  to 
become  greater  as  the  scour  hole  developed,  possibily  due 
to  a  strengthening  of  the  vortex  system  surrounding  the 
pile . 
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Fig.  Zf6.  -  Initial  stages  of  scour  devel 


opment 


?<4 


Fig.  W7*   -   Scour  pattern  30  minutes  after 
test  commencement 
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Fig.  Zf8.  -  Scour  pattern  2  to  3  hours  after 
test  commencement 
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Because  of  the  disturbance  of  flow  caused  by  the  pile, 
trailing  ripples  formed  immediately  behind  the  pile  ex- 
tending in  the  direction  of  wave  propagation  as  the  test 
progressed.   In  nearly  all  cases,  fairly  uniform  ripples 
formed  along  the  entire  sediment  bed.   The  wave  length  of 
these  ripples  ranged  from  approximately  1.5  to  2  inches. 
Most  of  the  ripples  that  covered  the  sediment  bed  initial- 
ly formed  at  the  interface  between  the  false  bottom  upstream 
of  the  pile,  and  the  sediment  bed.   The  ripples  probably 
formed  at  this  location  due  to  a  disturbance  in  flow,  and 
once  formed,  the  ripples  progressed  in  the  direction  of 
wave  propagation. 

In  all  cases,  the  maximum  scour  depth  occurred  down- 
stream of  the  pile,  approximately  1/4  inch  from  the  pile, 
and  at  an  angle  of  approximately  40°  to  the  left  or  right 
of  an  imaginary  line  drawn  through  the  pile  and  parallel 
with  the  wave  flume  walls. 
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PART  VII 
CONCLUSIONS 

1.  Scour  development  was  very  rapid  initially.   The  scour 
rate  then  decreased  until  an  equilibrium  scour  depth 
was  reached. 

2.  Scour  development  due  to  current  and  wave  action  oc- 
curred faster  than  scour  due  to  a  steady  current  alone. 

3 .  The  equilibrium  scour  depth  caused  by  wave-current 
interaction  is  approximately  10%  greater  than  the  scour 
depth  caused  by  a  steady  current  alone. 

4.  The  shape  of  the  scour  hole  due  to  waves  and  currents 
and  the  shape  of  the  hole  due  to  current  alone  are 
similar.   The  shape  resembles  an  inverted  cone. 

5.  For  a  constant  pile  diameter  and  sediment  size,  scour 
depth  increases  with  increasing  wave-current  velocity. 

6.  For  a  constant  wave-current  velocity  and  sediment  size, 
scour  depth  increases  with  an  increase  in  pile  diameter 

7.  For  a  current  velocity  greater  than  0.45  feet/ sec, 
scour  depth  increased  with  increasing  wave  height. 
For  a  current  velocity  less  than  0.45  feet/ sec,  scour 
depth  decreased  with  increasing  wave  height. 

8.  A  relationship  exists  between  scour  depth  and  sediment 
number  (N  ).   However,  scour  depth  is  not  dependent  on 
N  alone,  but  is  also  dependent  on  pile  diameter  and 
sediment  size. 


9.   A  relationship  also  exists  between  scour  depth  and 
the  pile  Reynolds  number  (NRp) .   Again,  scour  depth 
is  not  dependent  on  N^p  alone,  but  is  also  dependent 
on  pile  diameter,  as  well  as  sediment  size. 

10.  By  plotting  the  relative  scour  depth  (Sj/h)  against 
the  product  of  the  pile  Reynolds  number  and  sediment 
number,  (N^)  (N  )  ,  all  data  points  collapse  onto  a 
single  curve  regardless  of  variance  in  pile  diameter 
or  sediment  size. 

11.  A  relationship  exists  between  the  relative  scour  depth, 
(Sd/H1/3),  and  Njjp.Ng  and  (Ng)  (N^)  . 

12.  In  order  to  attain  the  same  scour  depth  with  a  single 
pile  as  that  achieved  in  a  multiple  pile  case,  the 
velocity  must  be  greater  assuming  all  other  variables 
are  identical. 

Although  the  amount  of  data  is  limited,  due  to  the 
fact  that  all  data  could  be  collapsed  onto  one  curve  when 
plotting  S,/h  versus  (Np-p)  (N  )  ,  (Figures  32  and  34)  ,  it 
appears  as  though  the  type  of  wave  train  used  for  model 
studies,  either  monochromatic  or  random,  is  insignificant 
in  the  prediction  of  scour  depth  about  a  single  cylindri- 
cal pile.   However,  when  attempting  to  predict  the  geometry 
of  the  scour  hole,  the  irregularity  of  the  wave  train  be- 
comes very  significant.   Lee  (13)  and  Armbrust  (2)  stated 
that  the  pattern  of  scour  for  wave-current  interaction 
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is  no  longer  conical  in  shape,  but  the  pattern  is  irregular 
due  to  the  wave  particle  orbital  motion.   The  data  in  this 
study  clearly  indicate  a  conical  scour  pattern  for  random 
wave-current  interaction.   Although  these  results  are  dif- 
ferent from  those  of  Lee  and  Armbrust,  they  agree  with 
field  data,  as  well  as  with  the  predictions  of  Niedoroda 
and  Dal ton  (15) .   The  fact  that  the  geometric  shape  of 
the  scour  hole  differs  with  the  type  of  wave  train  used, 
may  become  important  in  the  prediction  of  scour  depth 
around  multiple  pile  cases. 

Figures  (35)  and  (36)  show  the  relationship  found  for 
random  wave-current  interaction  between  the  relative  scour 
depth  Sj/H-,  i ^   and  the  sediment  number,  N  ,  and  pile  Rey- 
nolds number,  N„ p ,  respectively.   These  results  are  simi- 
lar to  those  found  by  Wells  and  Sorenson  (27)  and  present- 
ed in  Figures  (16)  and  (17) .   Because  of  similar  trends 
in  all  figures  and  due  to  the  fact  that  the  pile  Reynolds 

number  for  the  rapid  decrease  of  relative  scour  depth  is 

3 
approximately  7.5  x  10   in  both  cases,  (Figures  36  and  17), 

it  again  appears  as  though  the  irregularity  of  the  wave 
train  is  insignificant  in  the  prediction  of  scour  depth 
around  a  single  cylindrical  pile. 

In  addition,  Armbrust  (2)  found  that  there  was  an  in- 
crease of  scour  depth  with  increasing  wave  height.   As 
stated  earlier,  for  irregular  waves  when  the  current  ve- 
locity was  greater  than  .45  ft/ sec,  scour  depth  increased 
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with  increasing  wave  height,  and  the  scour  depth  decreased 
with  increasing  wave  height  for  a  current  velocity  less 
than  .45  ft/ sec. 

The  conclusions  resulting  from  this  study  were  found 
for  a  flow  in  which  the  steady  current  was  greater  than 
the  oscillatory  flow  produced  by  the  waves.   If  the  mag- 
nitude of  the  two  flows  in  relation  to  one  another  change, 
the  flow  patterns  will  also  change  and  the  geometry  of 
the  scour  hole  may  very  likely  differ  from  results  found 
in  this  study. 

As  mentioned  in  previous  investigations,  to  attempt 
to  predict  scour  depth  for  a  prototype  using  the  results 
of  this  study  could  lead  to  erroneous  results.  In  order 
to  accurately  predict  the  prototype  behavior,  complete 
similitude  is  required  between  model  and  prototype,  and 
field  data  is  required  to  evaluate  scale  effects,  should 
any  exist. 

Suggestions  for  Further  Research 

1.  This  study  could  be  repeated  under  the  same  conditions, 
thereby  omitting  scale  factors,  if  any  exist,  but  using 
monochromatic  waves  to  determine  if  the  irregularity 

of  the  wave  train  is  significant  in  scour  studies. 

2.  Further  research  can  be  conducted  to  investigate  the 
variation  of  scour  depth  with  varying  significant  wave 
height,  H, i-,  as  the  current  velocity  varies. 
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3.  It  would  be  desirable  to  investigate  why  there  is  a 

rapid  decrease  in  relative  scour,  S,/H-,  ,o,  for  a  pile 

3 
Reynolds  number,  N„p ,  of  approximately  7.5  x  10  . 

4.  The  study  could  be  repeated  using  various  sediment 

sizes  and  pile  diameters  to  determine  if  scale  effects 

are  negligible  when  using  a  plot  of  S,/h  versus  (N-np) 

(N  )  . 
v  sy 

5.  This  study  could  be  repeated  using  multiple  piles  in 
order  to  compare  the  results  with  previous  studies  in- 
vestigating scour  about  multiple  piles  due  to  monochro- 
matic wave-current  interaction. 


102 


REFERENCES 


1.  Abad,  G.N.,  and  Machemehl ,  J.L.,  "An  Experimental  Study 
of  Scour  Around  Marine  Foundations  Due  to  Oscillatory 
Waves  and  Unidirectional  Currents,"  Center  for  Marine 
and  Coastal  Studies  -  North  Carolina  State  University, 
Report  No.  74-5,  September,  1974. 

2.  Armbrust,  S.F.,  "Scour  About  A  Cylindrical  Pile  Due  to 
Steady  and  Oscillatory  Motion,"  Masters  Thesis,  Civil 
Engineering  Department,  Texas  A&M  University,  May,  1982. 

3.  Carstens,  M.R. ,  "Similarity  Laws  for  Localized  Scour," 
Journal  of  the  Hydraulics  Division,  ASCE,  Vol.  92,  No. 
HY3,  Proc.  Paper  4818,  May  1966. 

4.  Coleman,  N.L.,  "A  Theoretical  and  Experimental  Study  of 
Drag  Lift  Forces  Acting  On  A  Sphere  Resting  On  A  Hypo- 
thetical Stream  Bed,"  Proceedings,  Twelfth  Congress, 
IAHR,  Vol.  3,  Section  C18,  19,  1967. 

5.  Coleman,  N.L. ,  "Analyzing  Laboratory  Measurements  of 
Scour  at  Cylindrical  Piers  in  Sand  Beds,"  Proceedings , 
Fourteenth  Congress,  IAHR,  Vol.  3,  Section  C37,  1971. 

6.  Davies,  T.R.H.  and  Samad,  M.F.A.,  "Fluid  Dynamic  Lift 
On  A  Bed  Particle,"  Journal  of  the  Hydraulics  Division, 
ASCE,  Vol.  104,  No.  HY8 ,  August,  1978,  pp.  1171-1182. 

7.  Eagleson,  P.S.,  and  Dean,  R.G.,  "Wave  Induced  Motion  of 
Bottom  Sediment  Particles,"  Transactions ,  ASCE,  Vol.  126, 
1961,  Part  1,  pp.  1162-1189. 

8.  Einstein,  H.A.,  "Formulae  for  the  Transportation  of  Bed 
Load,"  Transactions,  ASCE,  Vol.  107,  1942,  pp.  561-577. 

9.  Graf,  W.H.,  Hydraulics  of  Sediment  Transport,  McGraw- 
Hill,  Inc.,  New  York,  1971,  pp.  83-86. 

10.  Hjorth,  P.,  "Studies  on  the  Nature  of  Local  Scour," 
Bulletin  Series  A  No.  46,  Department  of  Water  Resources 
Engineering  Lund  Institute  of  Techno logy/ University  of 
Lund,  Lund,  Sweden,  1975. 

11.  Jain,  S.C.  and  Fischer,  E.E.,  "Scour  Around  Bridge  Piers 
At  High  Flow  Velocities,"  Journal  of  the  Hydraulics 
Division,  ASCE,  Vol.  107,  No.  HY11,  November,  1980, 

pp.  1827-1840. 


103 


12.  Laursen,  E.M. ,  "Observations  on  the  Nature  of  Scour," 
Proceedings ,  Fifth  Hydraulic  Conference,  State  Uni- 
versity  of  Iowa,  Iowa  City,  Iowa,  1952,  pp.  179-197. 

13.  Lee,  T.Y. ,  "Scour  Around  A  Group  of  Circular  Piles 
Caused  By  Waves  and  Currents,"  Masters  Thesis,  Civil 
Engineering  Department,  Texas  A&M  University,  May, 
1983. 

14.  Newman,  J.N.,  Marine  Hydrodynamics,  The  MIT  Press, 
Cambridge,  Massachusetts,  1978,  pp.  237- 246. 

15.  Niedoroda,  A.W.  and  Dalton,  C,  "A  Review  of  the  Fluid 
Mechanics  of  Ocean  Scour,"  Ocean  Engineering,  Vol.  9, 
No.  2,  1982,  pp.  159-170. 

16.  Ninomiya,  K.  and  Tagaya ,  K. ,  "A  Study  on  Suction  and 
Scouring  of  Sit  On  Bottom  Type  Offshore  Structures," 
Proceedings,  Offshore  Technology  Conference,  OTC  1605, 
1972,  pp.  869-878. 

17.  Paintal,  A.S.,  "Concept  of  Critical  Shear  Stress  In 
Loose  Boundary  Open  Channels,"  Journal  of  Hydraulic 
Research,  IAHR,  Vol.  9,  No.  1,  1971,  pp.  91-109. 

18.  Palmer,  H.D. ,  "Wave-Induced  Scour  on  the  Sea  Floor," 
Proceedings ,  Civil  Engineering  in  the  Oceans  II,  ASCE, 
December,  1969,  pp.  703-716. 

19.  Ranee,  P.J.,  "The  Potential  For  Scour  Around  Large 
Objects,"  Seminar  on  Scour  Prevention  Techniques  Around 
Offshore  Structures ',    Society  of  Underwater  Technology, 
London,  England,  HT  December  1980. 

20.  Raudkivi,  A.  J.  and  Et tenia,  R.  ,  "Clear-Water  Scour  At 
Cylindrical  Piers,"  Journal  of  Hydraulic  Engineering, 
ASCE,  Vol.  109,  No.  3,  March,  1983,  pp.  338-350. 

21.  Shen,  H.W. ,  Schneider,  V.R. ,  andKaraki,  S.S.,  "Mechanics 
of  Local  Scour,"  Engineering  Research  Center,  Colorado 
State  University,  Fort  Collins,  Report  No.  CER66HWS22, 
June,  1966. 

22.  Shen,  H.W.,  Schneider,  V.R.,  and  Karaki,  S.S.,  "Local 
Scour  Around  Bridge  Piers,"  Journal  of  the  Hydraulics 
Division,  ASCE,  Vol.  95,  No.  HY6,  November,  1969,  pp. 
1919-1940. 

23.  Van  Hoften,  J.D.A.  and  Karaki,  J.,  "Interaction  of 
Wave  and  a  Turbulent  Current,"  Proceedings ,  Fifteenth 
Coastal  Engineering  Conference,  July  11-17,  1976, 
Honolulu,  Hawaii,  pp.  402-422. 


104 


24.  Vanoni,  V.  et  al.,  "Sediment  Transportation  Mechanics: 
Initiation  of  Motion,  Progress  Report  of  the  Task 
Committee  On  Preparation  of  Sedimentation  Manual,  Com- 
mittee on  Sedimentation,"  Journal  of  the  Hydraulics 
Division,  ASCE,  Vol.  92,  No.  HY2,  Proc.  Paper  4738, 
March  1966. 

25.  Venezian,  G.,  "Similitude  Relationships  and  Model  Test- 
ing," Dredging  Engineering  Shortcourse,  Paper  V-l, 
Texas  A&M  University,  January,  1984. 

26.  Wang,  R.K.  and  Herbich,  J.B.,  "Combined  Current  and 
Wave-Produced  Scour  Around  A  Single  Pile,"  Department 
of  Civil  Engineering  Texas  Engineering  Experiment 
Station,  COE  Report  No.  269,  Texas  A&M  University, 
March,  1983. 

27.  Wells,  D.R. ,  and  Sorensen,  R.M. ,  "Scour  Around  a 
Circular  Pile  Due  to  Oscillatory  Wave  Motion,"  Sea 
Grant  Publication  No.  208,  COE  Report  No.  113,  Texas 
A&M  University,  January,  1970 . 

28.  Yang,  C.T.,  "Incipient  Motion  and  Sediment  Transport," 
Journal  of  the  Hydraulics  Division,  ASCE,  Vol.  99, 

No.  HY10,  October,  1973. 


105 


APPENDIX  A 
Random  Wave  Spectra  Data 
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.6991 

.6516 

.6201 

.5978 

.5812 

.  5685 

.5587 

.5509 


TRANS2 


FINAL 


X  STROKE  =   129.659502619 


1.0343     9.9513 

1.0685     4.9788 

1.1028     3.2471 

1.1370     2.3317 

1.1713     1 

1 . 2055     1 

1 . 2398     1 

1.2741 

1 . 3083 

1 . 3426 

1 . 3768 

1.4111 

1.4453 

1.4796 

1.5138 

1.5431 
;804787138  mm 
mm  3.66471871869  STDEVS 


7572 
3764 
1327 
9907 
9146 
8747 
8538 
8435 
8400 
8412 
8458 
8523 


SIGNAL 

SPECTRUM 

mm-"-2/Hz 

0 .  00 

0 .  00 

1.57 

1643.48 

3629.65 

2176.26 

998. 10 

471. 13 

245.86 

140.47 

85.96 

55.49 

37.41 

26.  15 

13.85 

13.94 


I  SIGNAL  OUTPUT  =  6.51233980254  STDEVS 
■MING,  TRUNC  WAVES 
IGHTING  FUNCTION 

+0. 00000 
+25.71290 
-24.36528 
-4.92107 
+9 
+  1 


-1 


12731 
26716 
60015 
81417 
AT 


I  VALUE  STARTING 

5E   7F   65   39   14   F4   D7 
14   OF   08   03   00   FC   FS 

OCK  INTERVAL  =  .119047619048 


+42.63063 
+9. 11734 
-22. 24950 
+2.56339 
+9.09331 
-. 01803 
-2.79689 
-.65941 
ADDRESS  OOCO 
CA 
F3 
S  - 


+57.7463 
-5.5703 

-17.0950 
+7. 2982 
+7. 0203 


CF 
FA 
ID 


-1 


DA 
FA 
HEX 


7362 
7233 
6409 

E6 

FA 


F5 
FC 


FILTER 
FUNCTION 

0 . 0000 

0 . 0000 

.  0570 

1 . 8477 

2.7559 

2. 1434 

1.4592 

1 . 0087 

.7337 

.5589 

.4410 

.3577 

.2967 

.  2509 

.2156 

.  1378 


+45.7075 

-18.4574 

-11.7322 

+3.6747 


+3 


06 
FF 


6188 
50 1 9 
8087 
8160 

10 
FF 


MODEL 
FREQUENCY 
Hz 

.  1313 

.2625 

.3933 

.  5250 

.  656o 

.7375 

.9138 

1 . 0500 

1. 1313 

1.3125 

1.4433 

1.5750 

1 . 7063 

1 . 3375 

1.9688 

2. 1000 


13 

FE 


T I  MATED  GAIN  SETTING  =  7.33052823039 


EDBACK  OPTIONS: 

led back  Points  OA  03   OB  02   OF  01   11  03 

equence  Length   IK      2K      32K      12SK 


14 


03 

1M 


LOCITY  CONTROL 
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Y  = 


calculated  HS  =   9.98580007236 


^PLACEMENT  CONTROL 


WE 

WAVE 

TRANSFER  FUNCTIONS 

SIGNAL     FILTER     MODEL 

?EQUENCY 

SPECTRUM   TRANS 1 

TRANS2 

FINAL     SPECTRUM 

FUNCTION 

FREQUENCY 

Hz 

m-'-2/Hz 

mm- -2 /Hz 

Hz 

.01313 

0 . 0000 

9.6217 

1 . 0343 

9.9513 

0 .  00 

0. 

0000 

.  1313 

.  02625 

0 . 0000 

4.6596 

1 . 0635 

4.9788 

0 .  00 

0. 

0000 

.2625 

.  03938 

0 . 0000 

2.9445 

1. 1028 

3.2471 

0 .  00 

0. 

0000 

.3933 

. 05250 

0. 0000 

2.0507 

1. 1370 

2.3317 

0 .  00 

0. 

0000 

.  5250 

. 06563 

2.2426 

1 . 5002 

1. 1713 

1.7572 

69.25 

3807 

.6563 

. 07875 

32. 5300 

1. 1417 

1 . 2055 

1.3764 

616.28 

1. 

1406 

.7875 

. 09 1 88 

100.0229 

.9137 

1.2398 

1. 1327 

1283.41 

1. 

6547 

.9188 

.  10500 

125.7719 

.7776 

1.2741 

.  9907 

1234.49 

1. 

6328 

1 . 0500 

.  11813 

72.3776 

.6991 

1 . 3083 

.9146 

605.45 

1. 

1514 

1. 1313 

.  13125 

46.0842 

.6516 

1.3426 

.8747 

352.63 

8855 

1.3125 

.  14438 

31.3266 

.  620 1 

1 . 3768 

.8533 

228.35 

a 

7187 

1.4438 

.  15750 

22.2956 

.5978 

1.4111 

.8435 

153.62 

6047 

1 . 5750 

.  1 7063 

16.4187 

.5812 

1 . 4453 

.  8400 

115.85 

■ 

5221 

1 . 7063 

.  13375 

12.4146 

.5685 

1.4796 

.3412 

87.85 

■ 

4598 

1 . 8375 

.  19688 

9.5875 

.5587 

1.5138 

.8458 

68.  58 

■ 

4112 

1.9633 

.21 000 

7 .  5336 

.  5509 

1.5431 

.8528 

54.79 

3723 

2. 1000 

jDEV  PADDLE  DISPLMT  = 

25.2254078414  mm 

1  STROKE 

:  =   129.6595 

02619  mm  5.14003592362  STDEVS 

IX  SIGNAL 

.  OUTPUT  =  7. 

31375797316  STDEVS 

IIGHTING 

FUNCTION 

+0. 00000 

+42. 

31362 

+46 

.  0740 

+17.2287 

— 8 .  536  12 

-16. 

32912 

-16 

.5965 

-15.6599 

-8.43808 

+3. 

58199 

+  10 

.  0639 

+8. 19C 

>5 

+4.94569 

+3. 

34469 

+ 

.  1 1.29 

-3.9883 

-4. 13193 

-. 

9950 1 

+ 

.2543 

-.70S 

>6 

-.25819 

+  1. 

14111 

+ 

.3952 

-1. 19Z 

•o 

-. 35002 

+  1. 

47241 

+  1 

.  0045 

-.442 

a 

-.05661 

+  . 

71136 

- 

.  4055 

- 1 . 383 

■■? 

■:x  VALUE 

STARTING  AT 

ADDRESS 

OOCO 

>   76   7F 

'   2F   E9   D3 

D2   C 

>5   E9   OA 

1C 

1 7   OE 

09   00 

FZ 

3   FD   0 1 

FE   FF   03 

01   FD   FF   04 

03 

FF   00 

02   FF 

FC 

BCK  INTERVAL  =  .119047619048  S  =  ID  HEX 

3TIMATED 

GAIN  SETTING 

=  6.31005305647 

[EDBACK  OPTIONS: 

resdback 

Points  OA  03 

OB  02 

OF  01   11  0 

3   14  03 

sequence 

Length   IK 

2K 

32K 

12 

3K     1 M 

ILOCITY  CONTROL 
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APPENDIX  B 
Scour  Patterns  for  Random  Wave-Current  Interaction 


Ill 


.  ^I^.lf^^a.j-?f  ■-•■■■»;*■.  AkkAA*'^- 


Run  5. 


h 


=   1.167   ft 


Uc     =   .i+05   fps 

U       =    .Vf2    fps 

wave   spectrum  =  Darbyshire 

H1/3=    .2278   ft 


D 


=   .125    ft 


Sd     =    .126    ft 
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Run  6. 


h 


=  .917  ft 


UQ   =  ,330  fps 

U   =  .394  fps 

wave  spectrum  =  Darby shire 

H1  /,=  .1656  ft 

D 


=  .125  ft 

=  .115  ft 
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Run  7. 


h 


=   .917  ft 


Uc      =   .380  fpS 

U       =   .393   fps 

wave   spectrum  =   Pierson-Moskowitz 

H1/3=    .1266    ft 

D       =    .123    ft 

SH      =   .13   ft 


114 


Run  8. 


h       =  .917   fp 

Uc     =  .380  fps 

U       =  .401    fps 

wave   spectrum  =   Jonswap    (T  =   8  sec 

H1/3=  .1715   ft 

D        =  .125    ft 

Sd     =  .127   ft 


115 


Run  9. 


h       =  .917   ft 

Uc      =  .330   fps 

U       =  .386   fps. 

wave    spectrum  =   Jonswap    (T   =    10   s< 

H1/3=  .1930    ft 

D        =  .125    ft 

Sd     =  .09   ft 


116 


Run   10. 


h 


=   1.167   ft 


Uc     =   .^ 05   fps 

U       =    .k37   fps 

wave   spectrum  =   Pierson-Moskowitz 

H1/3=    .1317   ft 

D        =    .125    ft 

Sd      =    .122   ft 


117 


Run   1 1 . 


h 


=   1.167   ft 


Uc      =   .405    fps 

U       =   .459  fps 

wave   spectrum  =   Jonswap    (T  =   8  sec 

H1/3=    .1791    ft 

D       =   .175    ft 


SA      =    .152   ft 
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Run   12. 


h  =  1.167  ft 
Uc  =  .Zf05  fps 
U       =  .419  fps 

wave   spectrum  =   Jonswap    (T  =    10   se 

1/3: 


H,  /,=    .210    ft 


D 


=    .125    ft 


sd    =  .113  ft 
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Run  17. 


h   =  .917  ft 

Uc  =  ,650  fps 

U   =  .665  fps 

wave  spectrum  =  Darby shire 

H1/3=  .13^2  ft 


D 


=  .125  ft 


Sd  =  .18  ft 


Run   1  8. 


h        =    .917    ft 

Uc     =    .650   fps 

U°     =   .659  fps 

wave   spectrum  =   Pierson-Moskowitz 

H1/3=   .0967   ft 
D        =    .125    ft 
Sd     =    .174   ft 


121 


Run  19. 


h 


=  .917  ft 


Uc  =  .650  fps 

U   =  .660 

wave  spectrum  =  Jonswap  (T 

H1/3=  .1605  ft 

D   =  .125  ft 

Sd  =  .176  ft 


=  8  sei 
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Run  20. 


h       =   .917   ft 

Uc     =   .650   fps 

U       =   .   664   fps 

wave    spectrum  =  Jonswap(T  =    10   se 

H1/3=    .1093    ft 

D        =    .125    ft 


Sd      =    .168   ft 


APPENDIX  C 
List  of  Notations 
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Symbol 

a 
D 
D 


50 


1/3 


ES 


N. 
N 
N 


FP 
RP 


T 

U 

u 

( 

u 

] 
v 

P 


-  P 


(j0 

0 

T 

n 


Quantity 

wave  amplitude 

pile  diameter 

mean  diameter  of  sediment 

acceleration  due  to  gravity 

still  water  depth 

significant  wave  height 

wave  height 

wave  length 

scour  Euler  number 

flow  Froude  number 

pile  Froude  number 

pile  Reynolds  number 

sediment  number 

ultimate  scour  depth 

wave  period 

mean  wave-current  velocity 

current  velocity 

Ursell  parameter 

kinematic  viscosity  of  water 

fluid  density 

sediment/fluid  density 
difference 

angular  velocity  of 
oscillation 

angle  of  repose 

average  shear  stress 

water  surface  elevation 


Units   Dimensions 


ft 

L 

ft 

L 

ft 

L 

ft/sec 

L/T2 

ft 

L 

ft 

L 

ft 

L 

ft 

L 

ft 

L 

sec 

T 

ft/sec 

L/T 

ft/sec 

L/t 

— ■ 



ft2/s 

L2/T 

slug/ft2 

M/L3 

slug/ft2 

M/L3 

rad/sec 

1/T 

lb/ft2 

M/LT 

ft 

L 

/ 


I  ^ 


- 


Thesis 

E113 

c.l 


Eadie 

Scour  about  a  single 
cylindrical  pile  due  to 
random  waves  and  cur- 
rent . 


Thesis 
E113 

c.l 


Eadie 

Scour  about  a  single, 
cylindrical  pile  due  to 
random  waves  and  cur- 
rent . 


